AD-A037  192  ARMY  ENGINEER  WATERWAYS  EXPERIMENT  STATION  VICKSBURG  MISS  F/G  13/2 

LIQUEFACTION  POTENTIAL  OF  DAMS  AND  FOUNDATIONS.  REPORT  3.  DEVEL— ETC(U) 
FEB  77  G Y BALAOI*  8 ROHANI 

UNCLASSIFIED  WES-RR-S-76-2-3  NL 


1 1 
AO 

A037I92 

■ 

, ^ 

1 

, 

■ ^ 

- 

1 

i 

END 

DATE 

FILMED 

4-77 

RESEARCH  REPORT  S-76-2 


LIQUEFACTION  POTENTIAL  OF 
DAMS  AND  FOUNDATIONS 


DEVELOPMENT  OF  AN  ELASTIC-PLASTIC  CONSTITUTIVE 
RELATIONSHIP  FOR  SATURATED  SAND 


George  Y.  Baladi,  Behzad  Rohani 

Soils  and  Pavements  Laboratory 
U.  S.  Army  Engineer  Waterways  Experiment  Station 
P.  O.  Box  631/  Vicksburg/  Miss.  39180 

February  197^ 

Report  3 of  a Series 


PnMrad  for  Office/  Chief  of  Engineers/  U.  S.  Army 
Washington/  D.  C 20314 


Umisr  CWIS  31145 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PACE  fWhmn  Dmtm  Enimrtd) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 

Research  Report  S-76-2  ■ — ' 


LIQUEFACTION  TOTENTIAL  OF  J)AMS  AND JOUNDATIONSj 
"Report  3^  DEVELOPMENT  OF  JUI  JLASTI^-PLASTIC 
CONSTITUTIVE  RELATIONSHIP  FOR  SATURA'tED  SAND  • 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.  S.  Army  Engineer  Waterways  Experiment  Station 

Soils  and  Pavements  Laboratory 

P.  0.  Box  631,  Vicksburg,  Miss.  39l80 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  /TT 

Office,  Chief  of  Engineers,  U.  S.  Army  I/' 

Washington,  D.  C.  2031^* 


16.  DISTRIBUTION  STATEMENT  fo/ R*porO 


CT  OR  GRANT 


to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 


CWIS  311^4  5 


-IS.  SECURITY  CLASS.  (*0/ t/ii« 

Unclassified 


15*.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


19.  KEY  WORDS  (Contlnum  on  rmvmtmm  mldm  //  nmcmmamry  and  Idmntliy  by  block  numbor) 

Constitutive  relations 

Sands 

Dams 

Saturated  soils 

Foundations 

Liquefaction  (soils) 

ACT  fCantfoM  «a  fwprwm  If  nm—my  mad  fdmniHr  t»y  block  numibmr) 

This  report  describes  the  initial  attempt  at  developing  a three- 
dimensional  constitutive  relationship  for  saturated  sand.  An  isotropic  con- 
stitutive relationship  has  been  constructed,  within  the  framework  of  the  in- 
cremental theory  of  plasticity,  that  qualitatively  simulates  some  of  the 
features  of  the  stress-strain-pore  pressure  response  of  saturated  sand.  It 
accounts  for  the  hysteretic  behavior  of  pressure-volumetric  strain  response  of 

(Continued ) 


EDITION  OF  I NOV  6S  IS  OBSOLETE 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWlMn  Dmtm  BnltrtO) 


03 S'!  CO 


SECURITY  CLASSiriCATlON  OF  THIS  PAGE(irh«n  Oaf*  Bntmfd) 


20.  ABSTRACT  (Continued). 

^ sand,  the  effect  of  superimposed  hydrostatic  stress  on  shearing  response,  and 
the  shear-induced  volume  change.  It  does  not,  however,  treat  strain-softening 
behavior  and  does  not  predict  progressive  increase  of  pore  pressure  under  low- 
amplitude  (subyield)  cyclic  loading  conditions. 

The  behavior  of  the  constitutive  relationship  under  conventional  un- 
drained triaxial  test  condition  is  examined  and  correlated  with  experimental 
data  for  saturated  Reid  Bedford  Model  sand. 

A recommendation  is  made  that  the  present  constitutive  relationship  be 
extended  to  quantitatively  simulate  the  basic  features  of  the  stress-strain- 
pore  pressure  response  of  saturated  sand,  including  the  strain-softening  be- 
havior and  the  progressive  increase  of  pore  pressure  observed  under  cyclic 
loading  conditions. 


\ 


Unclassified 


security  classification  of  this  PAGEflW»*n  Dmtm  Bnfrmd) 


PREFACE 


This  investigation  was  conducted  by  the  U.  S.  Army  Engineer  Water- 
ways Experiment  Station  (WES)  under  the  sponsorship  of  the  Office,  Chief 
of  Engineers,  Department  of  the  Army,  as  a part  of  Project  CWIS  31ll*5, 
"Liquefaction  Potential  of  Dams  and  Foundations  Under  Earthquake 
Excitation. " 

The  investigation  was  conducted  by  Drs.  G.  Y.  Baladi  and  B.  Rohani 
during  the  period  September  1975  through  May  19T6,  under  the  general  di- 
rection of  Mr.  J.  P.  Sale,  Chief,  Soils  and  Pavements  Laboratory  (S&PL), 
and  Drs.  J.  G.  Jackson,  Jr.,  Chief,  Soil  Dynamics  Division  (SSD),  F.  G. 
McLean,  Chief,  Earthquake  Engineering  and  Vibrations  Division  (EE&VD), 
and  W.  F.  Marcuson,  III,  Project  Leader  (EE&VD).  The  report  was 
written  by  Drs.  Baladi  and  Rohani. 

COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE,  were  Directors  of 
WES  during  the  preparation  of  this  report.  Technical  Director  was 
Mr.  F.  R.  Brown. 
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DEVELOPMENT  OF  AN  ELASTIC-PLASTIC  CONSTITUTIVE 
RELATIONSHIP  FOR  SATURATED  SAND 

PART  I : INTRODUCTION 

Background 

1.  Liquefaction  may  be  defined  as  the  "transformation  of  a granular 

material  from  a solid  state  to  a liquefied  state. The  process  which 

leads  to  liquefaction  in  saturated  soil  is  believed  to  be  associated  with 

an  increase  of  pore  pressure  and,  therefore,  a concurrent  decrease  of 

intergranular  stresses.  Liquefaction,  and  subsequent  flow,  occur 

whenever  the  intergranular  stresses  and,  consequently,  the  shearing 

strength  of  the  material  are  substantially  reduced.  Experimentally  it 

has  been  shown  that  liquefaction  can  occur  under  both  monotonically 
2 3 U 

increasing  and  cyclic  ’ loading  conditions.  The  mechanism  of  liquefac- 
tion during  cyclic  loading,  however,  is  believed  to  be  somewhat  different 
than  that  under  monotonic-type  loading.^  For  example,  in  conventional 
triaxial  tests  (saturated  undrained  condition),  liquefaction  under  mono- 
tonic loading  can  occur  only  for  very  loose  sand.  Under  cyclic  loading, 
on  the  other  hand,  it  is  possible  to  achieve  liquefaction  for  dense  as 
well  as  loose  materials. 

2.  Regardless  of  the  loading  condition,  before  a rational  analysis 
of  a liquefaction  problem  in  saturated  sand  can  be  performed,  a constitu- 
tive relationship  describing  the  stress-strain-pore  pressure  response  of 
the  material  must  be  available.  The  constitutive  relationship  must  be 
expressed  in  three-dimensional  geometry  and  be  applicable  to  any  state 

of  stress  and  deformation.  Furthermore,  the  numerical  values  of  the 
parameters  (material  constants)  in  the  constitutive  relationship  should 
be  readily  derivable  from  meaningful  laboratory  test  data.  The  constitu- 
tive relationship  will  serve  two  purposes:  (a)  it  provides  a means  for 

the  interpretation  and  organization  of  laboratory  test  data  for  various 


states  of  stress  and  defonnation,  and  (b)  it  can  be  utilized,  in  con- 
junction with  the  field  equations,  to  perform  effective  stress  analysis 
for  dams,  or  other  earth  structures,  when  subjected  to  different  transient 
or  static-type  boundary  loading  conditions.  Once  a calculation  scheme 
for  effective  stress  analysis  has  been  established,  the  potential  for 
liquefaction  can  be  assessed  by  examining  the  stress-strain-pore  pressure 
distribution  within  the  earth  structure  of  interest.  It  should  be  pointed 
out  that  within  the  framework  of  this  type  of  analysis,  it  is  the  process 
that  leads  to  liquefaction  which  is  predicted  rather  than  the  actual  flow 
of  the  liquefied  material. 

3.  As  part  of  requirements  previously  mentioned,  a constitutive 
relationship  for  saturated  sand  must  be  able  to  treat  a two-phase  material 
consisting  of  solid  particles  and  pore  fluid.  Within  the  framework  of  the 
theory  of  continuous  mass  media,  which  is  the  basis  for  all  stress  and  de- 
formation analyses,  there  are,  in  principle,  two  approaches  which  can  be 
taken  to  formulate  a constitutive  equation  for  a two-phase  system.  The 
first  approach  is  based  on  the  continuum  theory  of  mixtures.^  The  second 

approach,  which  may  be  referred  to  as  a "pseudo  mixture  theory,"  is  based 
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on  the  elastic-plastic  theory  of  solids.  In  this  approach,  a two-phase 
system  is  simulated  by  simply  limiting  the  overall  compressibility  of  the 
material  to  that  of  the  pore  fluid  and  solid  particles  when  drainage  is 
not  allowed.  The  continum  theory  of  mixtures  in  the  first  approach  is 
quite  complicated  in  its  full  generality.  A simplified  version  of  the  mix- 
ture theory  has  been  applied  to  solution  of  liquefaction  problems  in  satu- 
0 

rated  sand.  However,  because  of  the  assumption  that  the  solid  skeleton  of 

thf  material  is  elastic,  the  theory  could  not  predict  progressive  increase 

9 

of  pore  pressure  under  cyclic  loading  conditions.  The  elastic-plastic 
constitutive  relations  are  generally  based  on  physical  observation  of 
material  behavior  under  laboratory  conditions,  and  have  been  used  success- 
fully for  the  solution  of  a number  of  geotechnical  problems They 
are  not  mathematically  too  constrained  and  one  could  easily  incorporate 
any  physically  reasonable  material  behavior,  at  least  in  an  ad  hoc  manner, 
in  their  formulation. 


Ob.lective 


U.  The  overall  objective  of  the  analytical  studies  of  the  lique- 
faction and  deformation  problem  of  saturated  sands  is  to  develop  a 
rational  calculation  framework  for  performing  effective  stress  analyses 
in  realistically  posed  boundary-value  problems.  As  was  pointed  out 
previously,  once  such  a calculational  framework  has  been  established,  the 
potential  for  liquefaction  under  a specified  boundary  loading  condition 
can  be  assessed.  The  investigation  is  to  be  carried  out  in  three  phases: 

a.  Phase  I.  The  first  phase  of  the  investigation  will  be 
devoted  to  formulation  of  a three-dimensional  elastic - 
plastic  isotropic  constitutive  model  for  saturated  sand. 
During  this  phase  attempts  will  be  made  to  incorporate 
in  the  model  some  of  the  basic  laboratory  observed  stress- 
strain-pore  pressure  behavior  of  saturated  sand.  The 
constitutive  model  will  then  be  examined  in  light  of  labo- 
ratory test  data  pertinent  to  the  liquefaction  phenome- 
non, at  least  qualitatively,  such  as  undrained  triaxial 
test  results,  and  any  shortcomings  of  the  model  will  be 
delineated. 

Phase  II.  During  the  second  phase  of  the  investigation  the 
constitutive  model  will  be  extended,  if  necessary,  so  that 
it  can  "successfully"  simulate  laboratory  behavior  of  satu- 
rated sand  pertinent  to  the  liquefaction  problem. 

c_.  Phase  III.  During  this  phase,  efforts  will  be  made  to  incor- 
porate the  constitutive  model  into  a suitable  numerical  com- 
puter code  for  subsequent  use  in  conducting  effective  stress 
analysis  (and  assessing  liquefaction  potential)  under  various 
boundary  loading  conditions. 

The  objective  of  this  report  is  to  document  the  results  of  Phase  I of  the 
investigation. 


Scope 

5.  The  development  of  the  constitutive  model  is  presented  in  Part  II. 
The  behavior  of  the  model  under  axisymmetric  triaxial  test  condition  is 
demonstrated  in  Part  III.  Part  IV  contains  the  results  of  comparison  of 
model  prediction  with  laboratory  behavior.  Conclusions  and  recommendations 
are  given  in  Part  V.  Appendix  A contains  the  fundamental  basis  of  elastic- 
plastic  constitutive  models  and  is  included  for  reference  purposes  and 
future  use. 
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PART  II:  DEVELOPMENT  OF  CONSTITUTIVE  MODEL 
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Mechanical  Behavior  of  Cohesionless  Soil 

6.  The  mechanical  response  of  a cohesionless  soil  when  subjected  to 

externally  applied  loads  is  a function  of  the  volumetric  and  deviatoric 

stress-strain  properties  of  the  material.  The  stress-strain  properties 

are,  in  turn,  affected  by  such  factors  as  void  ratio,  degree  of  saturation, 

interstitial  pore  fluid,  and  the  loading  history  of  the  material.  Void 

ratio,  in  particular,  reflects  the  state  of  compaction  of  the  material, 

e.g.,  loose  and  dense.  Loose  sand  compacts  and  exhibits  a ductile-type 

stress-strain  behavior  when  subjected  to  a deviatoric  state  of  stress. 

Dense  sand,  on  the  other  hand,  dilates  and  exhibits  a brittle-type  stress- 

strain  behavior  when  subjected  to  a similar  stress  condition.  The  boundary 

between  the  loose  and  dense  states  is  characterized  by  the  void  ratio  at 

which  shearing  deformation  occurs  without  volume  change.  This  void  ratio 

is  referred  to  as  "critical  void  ratio"  and  its  magnitude  varies  inversely 

7 

with  mean  normal  stress.  This  basic  difference  in  the  physical  behavior 
of  loose  and  dense  sand  seems  to  exist  only  at  relatively  low  confining 
pressures.  At  very  high  confining  pressures  (i.e.,  at  confining  pressures 
above  the  preconsolidation  pressure)  the  shearing  stress-strain  behavior 
of  dense  sand  also  resembles  a ductile-type  behavior  and  is  accompanied 
by  compaction.  Therefore,  the  stress-strain  behavior  of  sand  is  highly 
dependent  on  the  confining  pressure  or  the  superimposed  hydrostatic  state 
of  stress. 

7.  The  stress-strain  behavior  of  saturated  sand  also  varies  greatly 
depending  on  whether  sand  is  in  the  drained  or  undrained  state  during 
loading.  Since  the  pore  fluid  is  relatively  incompressible,  in  the  un- 
drained condition  pressure  builds  up  in  the  pore  fluid.  The  effective 
stress  carried  by  the  soil  particles  is  the  total  stress  minus  the  pore 
pressure.  It  is,  therefore,  the  effective  confining  pressure  that 
influences  the  stress-strain  behavior  of  saturated  sand. 

8.  In  summary,  the  stress-strain  behavior  of  a satiirated  sand 
depends  on  its  initial  void  ratio,  the  current  and  past  loading  history. 


T 


1 


and  the  effective  confining  pressure.  Yield  and  deformation  in  the 
undrained  condition  typically  resemble  the  characteristics  shown  in 
Figure  1.  In  the  next  section,  the  development  of  an  elastic-plastic 
constitutive  relationship  that  partially  describes  the  behavior  shown  in 
Figure  1 is  presented. 

Constitutive  Model 

9.  The  fundamental  basis  of  elastic-plastic  constitutive  models 
is  presented  in  Appendix  A.  The  elastic  behavior  of  the  models  is  de- 
fined by  Equations  All  through  Al6  of  Appendix  A.  The  plastic  behavior 
of  the  models  is  described  by  Equations  AIT  through  A28 . The  complete 
elastic-plastic  description  is  governed  by  Equations  A29  and  A30.  To 
apply  these  equations  to  a two-phase  continuum,  consisting  of  solid 
skeleton  and  pore  fluid,  the  normal  stress  components  should  be  divided 
into  two  parts;  the  stress  carried  by  the  solid  structure,  referred  to  as 
effective  stress,  and  the  stress  carried  by  pore  fluid,  referred  to  as 
pore  water  pressure.  Mathematically,  total  stress  can  be  expressed  as 


•.  I 


0 = 0 + u 


(la) 


or,  in  tensorial  form,  as 


f ' 


it 


a. . = a. . + u 6.  . 

aj  11  11 


(lb) 


where 


'^iJ 


ij 


= total  stress  tensor 
= effective  stress  tensor 
= Kronecker  delta 


For  example,  in  the  case  of  triaxial  test  in  cylindrical  coordinate  system 
z , r , and  0 , Equation  1 takes  the  following  form 


a + u 
z 


+ u = Cq  + u 


(2) 


where  and  are,  respectively,  the  axial  and  radial  stress 

components.  The  mathematical  forms  of  various  response  functions  for  the 
proposed  constitutive  relationship  for  saturated  sand  are  presented  in 
the  following  paragraphs . 

10.  For  the  elastic  (recoverable)  response  of  the  material  it  is 
assumed  that  the  bulk  modulus  K is  a linear  function  of  pressure;  thus. 


K . Kj  + j’/3  - K.  t 


(3) 


where  K.  and  K,  are  material  constants  which  can  be  determined  from 
1 1 

the  slope  of  the  unloading  curve  from  an  isotropic  consolidation  test 

f 

(Figure  2).  The  q^uantity  is  the  first  invariant  of  stress  tensor 

» r 


and  P is  the  effective  mean  normal  stress  (P 


for  triaxial 


test).  Tlie  elastic  shear  modulus  G is  assumed  to  be  constant  (in  general 
the  shear  modulus  is  a function  of  the  state  of  stress  and  strain.  For 
simplicity  and  as  a first  order  of  approximation,  however,  it  is  assumed 
to  be  constant ) . 

11.  For  the  plastic  behavior  the  loading  function  (Equation  A9) 
is  assumed  to  be  isotropic  and  to  consist  of  two  parts  (Figure  3):  an 

ultimate  failure  envelope  which  serves  to  limit  the  maximum  shear  strength 
of  the  material 


and  a strain-hardening  cap 


F , k)  = (J^  - hf  + - (X  - hf 


(U) 


(5) 


The  quantity  is  the  second  invariant  of  stress  deviation  tensor 


(j„  = (a  - o ) /3  for  triaxial  test).  The  hardening  function  k which 
£:  z r 

controls  material  compaction  and/or  dilatation  is  assumed  to  be  equal  to  the 

p 

plastic  volinnetric  strain  ia-kes  the  following  form 
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(6) 


P 1 , ^2  X ^ . 1 . X ^ , s 

\k  Kg  Kq  3 ” K^  3 ^ 

t 

The  intersection  of  the  cap  with  the  axis  is  denoted  hy  X . The 

ratio  of  the  major  to  the  minor  axis  of  the  elliptic  cap  is  denoted  by 

T 

R . The  value  of  at  the  center  of  the  cap  is  denoted  by  L which 

is  related  to  other  parameters  by  the  relation 


Equation  7 indicates  that  the  failure  envelope  (Equation  U)  intersects 
each  ellipse  at  the  crown. 

12.  The  material  constants  and  in  Equation  6 are  deter- 

mined from  the  slope  of  the  loading  curve  from  an  isotropic  consolidation 
test  (Figure  2).  The  conditions  of  uniqueness  and  stability  have  to  be 
satisfied  in  determining  the  material  constants  ’ ^0  ’ 

K„  . These  conditions  are  satisfied  if  the  following  inequalities  are 


adhered  to: 


^i 


KilKg 


The  parameter  M (Equation  h)  is  indicative  of  the  frictional  strength  of 
the  material  and  is  related  to  the  effective  angle  of  internal  friction 
ij)  through  the  following  relation  (for  triaxial  compression) 


M = (9) 

3 - sin  <() 

13.  In  summary,  there  are  seven  material  constants  that  describe 
the  behavior  of  the  proposed  model,  i.e.,  , Kg  , G , M , 

and  R . A more  complicated  (and  possibly  more  realistic)  model  can  be 
developed  by  replacing  these  constants  with  appropriate  mathematical 
expressions  that  are  functions  of  stress  and/or  strain  invariants. 

1^4 . It  should  be  noted  that  the  proposed  model  allows  for  the  el- 
liptic caps  to  expand  as  well  as  to  translate  relative  to  the  origin  of 


the  .1^  ’ axis.  Thus,  the  model  incorporates  both  isotropic  and 

some  form  of  kinematic  hardening.  In  the  special  case  where  X = 2L  , 


i.e.,  R = , the  strain-hardening  cap  is  not  allowed  to  translate 

and  only  expands  isotropically.  In  this  case,  the  number  of  material 
const  .its  reduces  to  six  and  the  plastic  behavior  of  the  model  becomes 
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similar  to  the  model  proposed  by  Schofield  and  Wroth.  In  the  present 
model,  however,  only  a portion  of  the  elliptic  cap  is  used  as  loading 
function,  i.e.. 


(^^2  » '^1  ^ . 


if  L ^ 


■ 4 • 4>  If 


Schofield  and  Wroth,  on  the  other  hand,  use  the  entire  cap  as  loading 

function.  It  is  believed  that  the  use  of  the  entire  cap  violates  the 
12 

Drucker  stability  postulate  which  is  sufficient,  although  not  necessary, 
to  satisfy  all  thermodynamic  and  continuity  requirements  of  the  incremental 
theory  of  plasticity.  Stability  ensures  that  all  physically  reasonable 
initial -boundary-value  problems  are  properly  posed  in  the  mathematical 
sense . 

15.  An  undrained  condition  is  simulated  by  the  model  by  imposing 
the  condition 

-0  (11) 

where  is  ihe  increment  of  total  volumetric  strain.  The  model  then 

calculates  the  stress  path  (and  the  associated  material  response) 
corresponding  to  Equation  11.  This  stress  path  is  assumed  to  be  the 
effective  stress  path  that  the  material  will  experience  during  an  un- 
drained test.  Within  the  framework  of  elastic-plastic  models  (Appendix  A) 
Equation  11  can  be  satisfied  by  allowing  the  increment  of  plastic  volumetric 

p 

strain  “to  be  the  negative  of  the  increment  of  elastic  volumetric 


the  response  of  the  material  (including  volumetric  strain)  for  any  i 

specified  stress  path.  The  response  of  the  material  under  axisymmetric  i 

) 

triaxial  test  condition  is  demonstrated  in  Part  III.  \ 
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PART  III : BEHAVIOR  OF  THE  CONSTITUTIVE  MODEL  IN 

AN  AXISYMMETRIC  TRIAXIAL  TEST  CONDITION 


l6.  Most  of  the  mechajiical  testing  of  sand  for  engineering  purposes 
is  performed  in  the  triaxial  test  apparatus.  It  is  of  interest,  therefore, 
to  investigate  the  behavior  of  the  proposed  model  in  a triaxial  test  under 
both  drained  and  undrained  conditions.  Adopting  the  z-axis  of  a cylindri- 
cal coordinate  system  (z  , r , and  0)  as  the  axis  of  symmetry  of  the  soil 
sample,  the  effective  stress  tensor  and  strain  tensor  associated  with  the 

I I 

triaxial  test  become  (a.  = a and  = e ) 

0 r 0 r 


o'  0 0 

z 


(12a) 


e . = 0 


(12b) 


The  variables  P = Jj/3  (effective  mean  normal  stress),  (the  second 

invariant  of  stress  deviation  tensor),  and  (mean  volumetric  strain) 

KK 

associated  with  the  above  stress  and  strain  tensors  take  the  following 
forms 


P = J^/3 


o'  + 2a' 

z r 


(13a) 


■>2 -I 


(13b) 


e + 2e 
z r 


(13c) 


i. 


;l 


Isotropic  Consolidation  Phase 


17*  During  isotropic  consolidation  phase  of  the  triaxial  test 


o'  = o'  = J /3  = P 


(ll4a) 


kk 

e = e = —r— 
z r 3 


ilkh) 


The  relation  hetween  the  elastic  volumetric  strain  increment  and  the 
increment  of  effective  mean  normal  stress  is  given  as  (see  Equation  A12) 


dP  = K de. 


kk 


(15) 


where  the  elastic  hulk  modulus  K is  given  in  Equation  3*  Substituting 

Equation  3 in  Equation  15  and  integrating  the  resulting  expression  provides 

E 

the  following  relation  between  the  elastic  volumetric  strain  e and 

K.K 

pressure 


I" 


'' ' ! 

r I 

K • ■ 

if 

> I *. 

. 

r ,*  ‘ 

: 1 


El  f^l  ' 

"kk  = iT  K~  ^ 1 

1 '•1 


(16) 


The  relation  between  the  plastic  volumetric  strain  and  pressure  is 

iCK 

given  in  Equation  6.  During  isotropic  consolidation  — in  Equation  6 

f 3 

should  be  replaced  by  P , thus 


"kk  Kg 


In  view  of  Equations  l6  and  17  the  total  volumetric  strain  takes  the 
following  form 


jf 


1 (^2  . 'I 


) 


(18) 


Equations  l6  through  l8  provide  a complete  specification  for  the  behavior 
of  the  material  during  isotropic  consolidation  test. 
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l8.  The  qualitative  behavior  of  the  model  during  isotropic  consoli- 
dation is  shown  in  Figure  U.  The  slope  of  the  pressure-volumetric  strain 
curve  during  virgin  loading  can  be  obtained  from  Equation  l8 


= K2P  +k^  = k 


Combining  Equations  19  and  17  results  in 


K P + K. 

1 1 


1 + (K  P + K ) 

^ dP 


1 + K 


The  second  term  in  the  denominator  of  Equation  20  produces  an  apparent 

'X, 

softening  of  the  modulus  K due  to  plastic  compaction.  For  a very 
dense  sand,  the  softening  term  is  very  small,  and  the  modulus  K 
approaches  the  elastic  bulk  modulus  K . Also,  if  a sample  is  consolidated 
from  point  1 to  point  2 (Figure  U),  unloaded  from  point  2 to  point  3,  and 
then  reloaded  from  point  3 to  point  2,  the  model  dictates  that  the  reload- 
ing behavior  is  purely  elastic  and  the  second  term  in  the  denominator  of 
Equation  20  is  zero.  This  type  of  behavior,  however,  may  not  be  completely 
true  for  an  actual  sand  and  is  only  a mathematical  idealization. 

Shear  Phase 

19.  During  the  shear  phase  of  a conventional  triaxial  test  cell 
pressure  is  constant 


a = constant  = P 
r c 


(21a) 


da  = 0 
r 


(21b) 


where  P^  is  the  confining  pressure  at  the  end  of  the  isotropic  con- 
solidation phase.  For  an  undrained  test  the  volumetric  strain  is  constant 
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during  the  shear  phase  and  is  equal  to  the  volumetric  strain  achieved  at 
the  end  of  the  consolidation  phase  of  the  test.  Accordingly,  for  an 
undrained  test 


V 

A 


(22a) 


which  leads  to 


de  * - 77  de 
r 2 z 


(22b) 


1 

>■  't 

^ i 


The  plastic  volumetric  strain,  or  the  hardening  function,  during  the  shear 
phase  is  given  by  Equation  6.  In  view  of  Equations  l6  and  6,  the  total 
volximetric  strain  becomes 


o' 

1 1 „ 

X 

1 

+ 1 

+ — «,n 

1 Kg 

O |l 

1 ' 


-r- 

1 '•1 


(23) 


. < 
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For  the  \indrained  condition,  since  e is  constant  and  is  equal  to  the 

iCK 

volumetric  strain  achieved  at  the  end  of  the  isotropic  consolidation  phase 
(Equation  l8).  Equation  23  becomes 


'kk 


^2  1 

— P + li  = constant 
Kq  c ) 


P = P 


+ 1 


- (ic-  3 ^ ^ 


(2M 


The  equation  of  the  elliptic  cap  (Equation  5)  for  triaxial  configuration 
takes  the  following  form 


' L 2 r2  2 2 


(25) 
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Equations  2h  and  25  can  be  combined  to  determine  an  expression  for  the 

I 

effective  stress  path  (i.e.,  versus  P ) for  the  undrained  con- 

dition. Having  determined  the  effective  stress  path  the  pore  water  pres- 
sure can  be  obtained  from  Equation  1.  The  total  strain  can  be  obtained 
from  Equation  A29.  A computer  program,  called  TDRIVER,  was  developed  to 
numerically  solve  the  above  system  of  equations  and  generate  various  plots 
of  stress-strain-pore  pressure  response  for  undrained  conditions.  This 
program  and  its  flow  charts  are  available  upon  request.  Typical  results 
from  computer  code  TDRIVER  are  presented  in  Part  IV.  Figure  5 depicts 
qualitatively  the  effect  of  the  parameter  R on  the  stress-strain-pore 
pressure  response  for  conventional  triaxial  test.  It  is  noted  from 
Figure  5 that  R has  a substantial  influence  on  the  behavior  of  the 
model  (the  value  of  R reflects  the  relative  density  of  the  material, 
i.e.,  the  denser  the  material  the  smaller  the  value  of  R). 

20.  In  the  case  of  drained  test  the  response  of  the  material,  for 
a given  stress  path,  is  determined  from  Equations  23,  25,  and  A29.  Typical 
results  for  drained  condition  predicted  by  the  model  are  shown  qualitatively 
in  Figure  6. 


PART  IV:  COPfPARISON  OF  LABORATORY  TEST  DATA  WITH  MODEL  PREDICTION 

Experimental  Program 

21.  The  experimental  program  consisted  of  two  series  of  tests  con- 
ducted on  saturated  samples  of  Reid  Bedford  Model  sand.  The  gradation 
curve  for  the  Reid  Bedford  Model  sand  is  shown  in  Figure  7-  The  material 
consists  of  subround  to  subangular  particles.  Each  test  series  consisted 
of  a load/unload  isotropic  consolidation  test  and  four  conventional  con- 
solidated undrained  triaxial  tests.  Saturation  of  the  specimens  was 
achieved  by  backpressure  saturation.  The  first  test  series  was  conducted 
on  specimens  having  a relative  density  of  approximately  76  percent.  Rela- 
tive density  for  the  second  series  of  tests  was  about  38  percent.  Test 
data  from  the  first  series  of  tests  are  documented  in  Plates  1 through  3* 
Plate  1 depicts  the  results  of  the  isotropic  consolidation  test  presented 
in  terms  of  a plot  of  effective  mean  normal  stress  versus  volumetric  strain. 
The  results  of  the  triaxial  tests  are  shown  in  Plates  2 and  3 in  terms  of 
principal  stress  difference  and  excess  pore  pressure  versus  axial  strain, 
respectively.  The  corresponding  results  from  the  second  test  series  are 
documented  in  Plates  4 through  6. 

Material  Constants 

22.  As  was  pointed  out  previously,  there  are  seven  material  con- 
stants associated  with  the  proposed  constitutive  model  which  must  be  deter- 
mined experimentally.  Four  of  the  material  constants  (K^  ’ ^1  ’ *^0  ’ 

Kg)  are  associated  with  the  isotropic  consolidation  test  (see  Figure  2). 

The  parameter  M is  the  slope  of  the  ultimate  failure  envelope  of  the 
material  plotted  in  the  principal  stress  difference-effective  mean  normal 
stress  space  (see  Figure  3).  The  elastic  shear  modulus  G (assumed  to 
be  a constant  in  the  proposed  model)  can  be  estimated  from  unloading 
slopes  of  principal  stress  difference-principal  strain  difference  plots. 

The  parameter  R (the  ratio  of  the  major  to  the  minor  axis  of  the  ellip- 
tic cap)  can  be  determined  directly  from  test  results.  For  this 

study,  however,  since  K test  results  were  not  available,  the  value  of 


18 


the  parameter  R was  selected  in  order  to  obtain  a good  fit  to  the  experi- 
mental effective  stress  paths.  The  value  of  shear  modulus  was  selected  to 
obtain  a reasonable  fit  to  experimental  stress-strain  curves.  The  final 
values  of  the  material  constants  for  the  two  test  series  are  given  in 
Table  1. 


Comparison  of  Test  Results  with  Model  Behavior 

23.  Figures  8 through  10  depict  comparison  plots  of  test  results  -luth 
model  behavior  for  the  first  test  series  (relative  density  76  percent)  for 
consolidation  pressures  of  0.72  and  7*2  tsf.  Figure  8 shows  plots  of  effec- 
tive stress  paths  in  the  principal  stress  difference-effective  mean  normal 
stress  space.  Principal  stress  difference  versus  axial  strain  relations 
are  shown  in  Figure  9«  Figure  10  depicts  plots  of  excess  pore  pressure 
versus  axial  strain  relations.  The  corresponding  set  of  plots  for  the  sec- 
ond test  series  (relative  density  38  percent)  is  presented  in  Figures  11 
through  13.  For  the  purpose  of  model  calculations  the  yield  function  F 
(the  elliptic  cap)  was  used  as  loading  function  from  the  beginning  of  the 
test  until  the  effective  stress  path  reached  the  ultimate  failure  envelope 
f . From  then  on,  the  ultimate  failure  envelope  was  used  as  loading  func- 
tion. During  the  first  part  of  the  calculation  (i.e.,  using  F as  loading 
function)  excess  pore  pressure  increased  and  reached  its  maximum  value  when 
the  effective  stress  path  reached  the  ultimate  failure  envelope.  From  then 
on,  excess  pore  pressure  decreased  and  eventually  became  negative  as  the 
test  continued.  As  observed  from  Figures  8 through  I3,  the  proposed  consti- 
tutive model  qualitatively  simulates  the  stress-strain-pore  pressure  response 
of  the  material  for  both  test  series. 

2k.  In  order  to  evaluate  the  effect  of  the  parameter  R on  the 
stress-strain-pore  pressure  response  of  the  material  and  to  show  the  capa- 
bility of  the  model  in  simulating  the  behavior  of  loose  sand,  three  calcu- 
lations were  performed  for  a hypothetical  material  having  properties 

= TO  tsf  , = 37  , Kq  = ItO  tsf  , = 3.7  , M = 1.02  , and  G = 50  tsf. 

The  results  of  the  calculations  are  presented  in  Figure  llj  for  three 
different  values  of  R . The  curves  marked  B are  associated  with  a 
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base  line  value  of  R = = 5-09  • This  value  of  R corresponds  to 

D 

3*^ 

— which,  as  was  pointed  out  previously,  allows  for  the  strain-hardening 
cap  to  expand  isotropically  only.  The  curves  marked  A and  C are  associated 
with  R values  of  1+R  and  R„/^  , respectively.  It  is  noted  from 

D D 

Figure  1I4  that  the  parameter  R significantly  controls  the  stress-strain- 
pore  pressure  response  of  the  material.*  In  general,  for  a given  confin- 
ing pressure,  the  greater  the  value  of  R the  smaller  the  shear  strength 
of  the  material  and  the  larger  the  induced  pore  pressure.  For  a given  con- 
fining pressure,  smaller  values  of  R are  associated  with  higher  shear 
strengths  and  lower  pore  pressures . 


* For  the  hypothetical  material  the  calculations  were  terminated  when 
the  effective  stress  path  reached  the  ultimate  failure  envelope  f , 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


25.  A three-dimensional  elastic-plastic  isotropic  constitutive 
relationship  has  been  developed  that  quaMtatively  simulates  certain 
characteristics  of  the  stress-strain-pore  pressure  response  of  saturated 
granular  materials.  In  particular,  the  constitutive  relationship  accounts 
for  the  hysteretic  behavior  of  pressure-volumetric  strain  response  of 
sand,  the  effect  of  superimposed  hydrostatic  stress  on  shearing  response, 
and  the  shear-induced  volume  change. 

26.  The  constitutive  equation  does  not  treat  work-softening  behavior 
(Figure  1,  curve  No.  2),  and  can  not  predict  a progressive  increase  of 
pore  pressure  under  low-amplitude  (subyield)  cyclic  loading  conditions. 

Both  of  these  phenomena  are  important  features  which  contribute  to  lique- 
faction of  saturated  sand. 

27.  It  is  recommended  that  during  the  second  phase  of  this  investi- 
gation the  present  constitutive  equation  be  extended  to  account  for  the 
strain-softening  behavior  and* the  progressive  increase  of  pore  pressure 
observed  under  cyclic  loading  conditions.  This  extension  is  a major  task 
from  a theoretical  point  of  view,  but  can  be  accomplished  within  the  basic 
structure  of  the  present  model.  Preliminary  investigations  have  indicated 
that  such  an  extension  may  be  accomplished  in  a number  of  ways.  It  may  be 
necessary  to  try  several  techniques  and  adopt  the  one  which  best  simulates 
the  experimental  data  while  at  the  same  time  satisfying  all  theoretical 
requirements.  However,  as  a basic  requirement  to  any  feasible  technique 
the  compressibility  of  pore  fluid  and  soil  skeleton  must  be  included  in 
the  analysis  of  undrained  behavior. 

28.  If  the  extension  of  the  present  model  is  successfully  accom- 
plished, it  is  further  recommended  that  during  Phase  II  of  this  study  the 
shear  modulus  G and  the  parameter  R (the  ratio  of  the  major  to  the 
minor  axis  of  the  eliiptic  loading  function)  be  replaced  by  appropriate 
functional  forms  in  order  to  quantitatively,  as  well  as  qualitatively, 
simulate  the  stress-strain-pore  pressure  response  of  saturated  granular 
materials.  This  modification  is  a relatively  simple  task  (in  comparison 


rath  extension  of  the  model  discussed  in  paragraph  27)  from  a 
point  of  viev;.  The  particular  functional  forms  of  G and  R 
must  he  selected  based  on  experimental  observation  of  material 
under  diverse  loading  conditions. 
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Table  1 

Numerical  Values  of  Material  Constants 
for  Reid  Bedford  Model  Sand 


[ 

Test  Series 

K. 

1 

tsf 

"l 

"o 

tsf 

"2 

M 

G» 

tsf 

R* 

i 

i 

1 

1 

1 

(relative  density 

90.0 

2UU.0 

72.0 

160.0 

1.53 

2IO.O 

0.05 

= 16%) 

! ' 2 

(relative  density 

57.5 

208.0 

36.0 

IU2.O 

1.33 

100.0 

1.5 

= 38!?) 


* Selected  to  give  good  fit  to  experimental  data,  in  lieu  of  being 
determined  from  separate  tests. 
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Figure  2.  Proposed  relationship  for  isotropic  consolidation  test 


conventional  triaxial  test 


Drained  behavior  of  the  model  in  conventional  triaxial  and  constant 
mean  normal  stress  tests 
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Figure  8.  Comparison  of  effective  stress  paths  produced  by  the  model 
vri.th  test  data  for  saturated  Reid  Bedford  Model  sand  in  undrained 
triaxial  shear  (Test  Series  l) 
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Figure  9>  Comparison  of  principal  stress  difference  vs  axial  strain 
relations  produced  by  the  model  with  test  data  for  saturated  Reid 
Bedford  Model  sand  in  undrained  triaxial  shear  (Test  Series  l) 


Figure  10.  Comparison  of  excess  pore  water  pressure  vs  axial  strain 
produced  by  the  model  with  test  data  for  Reid  Bedford  Model  sand  in 
undrained  triaxial  shear  (Test  Series  l) 
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Figure  12.  Comparison  of  principal  stress  difference  vs  axial  strain 
relations  produced  by  the  model  with  test  data  for  saturated  Reid 
Bedford  Model  sand  in  undrained  triaxial  shear  (Test  Series  2) 
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Figure  13 . Comparison  of  excess  pore  water  pressure  vs  axial  strain 
produced  by  the  model  with  test  data  for  Reid  Bedford  Model  sand  in 
undrained  triaxial  shear  (Test  Series  2) 
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Figure  lU.  The  effect  of  the  parameter  R on  undrained  behavior  of  the  model  in  conventional 

triaxial  test  for  a hypothetical  material 
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APPENDIX  A 


FUNDAMENTAL  BASIS  OF  ELASTIC-PLASTIC  MATERIAL  MODELS 
Basic  Concepts  From  Continuum  Mechanics 

1.  In  engineering  practice  it  is  convenient,  and  often  reasonable, 
to  disregard  the  structural  details  of  materials  and  consider  their  gross 
behavior  only.  Engineering  materials  are  therefore  described,  or  charac- 
terized, mathematically  within  the  frameworks  of  the  theory  of  continuous 
mass  media.  Neglecting  thermal  effects,  the  basic  field  equations  that 
govern  the  motion  of  a continuum  are  the  continuity  equation* 

If-  + (p  V.  ) . = 0 (Al) 

do  1 


and  the  equations  of  motion 


+ F.  - 

1 


pa^  = 0 


(A2) 


where 

p = mass  density 
t = time 

V.  = components  of  velocity  vector 

a.  = symmetrical  stress  tensor 
J 

F^  = components  of  body  force 

a^  = components  of  acceleration  vector 

2.  Equations  Al  and  A2  constitute  four  equations  that  involve  ten 

unknown  functions  of  time  and  space:  the  mass  density  p , the  three 

velocity  components  v.  , and  the  six  independent  stress  components  a. 

r 1 j 

* Indices  take  on  values  1,  2,  or  3.  A repeated  index  is  to  be  summed 
over  its  range.  Comma  in  the  subscripts  represents  a derivative. 
Quantities  are  referred  to  rectangular  Cartesian  coordinates  X.  . 


Al 


The  body  force  components  are  known  quantities  and  the  acceleration 

components  a^  are  expressible  in  terms  of  the  velocity  components  . 

Therefore,  six  additional  equations  relating  the  ten  unknown  variables 
are  required  in  order  to  determine  the  motion  or  deformation  of  a medi- 
um when  subjected  to  external  disturbances  such  as  surface  forces  and/or 
displacements.  In  continuum  mechanics,  such  relations  are  stated  by  con- 
stitutive equations  (or  material  models),  which  relate  stresses  to  defor- 
mation and  history  of  deformation.  The  difference  between  constitutive 
equations  and  field  equations  (Equations  Al  and  A2)  is  that  the  latter 
are  applicable  to  all  materials,  whereas  the  former  represent  the  intrinsic 
response  of  a particular  material  or  class  of  materials. 

3.  The  general  form  of  a constitutive  equation  may  be  expressed 
by  the  functional  form 


1 > P)  = 0 (A3) 

ab  mn  rs  qp  ij 

where  the  deformation-rate  and  spin  tensors,  D and  , respectively, 

mn  qp 

are  related  to  the  components  of  the  velocity  vector  v^ 


D = i (v  + V ) 
mn  d m,n  n,m 


£1  = — (v  - V ) 

pq  2 p,q  q,p 


and  the  strain  tensor  e is  related  to  the  components  of  displacement 

rs 

vector  Uj^  . For  small  displacement  gradients 


e = r (u  + u ) 
rs  2 r,s  s,r 


Equations  Al  through  A3  constitute  ten  equations  which  include  ten  unknown 
variables.  These  equations  will  lead,  in  conjunction  with  the  kinematic 
relations  given  by  Equations  kk  and  A5  and  boundary  conditions , to  a com- 


plete description  for  solution  of  a boundary-value  problem 


h.  The  mechanical  behavior  of  a number  of  engineering  materials 
is  described  within  the  framework  of  elastic-plastic  constitutive  relation- 
ships. The  development  of  the  specific  functional  form  of  Equation  A3 
for  the  elastic-plastic  models  is  given  in  the  following  section. 

General  Description  of  Elastic-Plastic  Constitutive  Models 

5.  The  basic  premise  of  elastic-plastic  constitutive  models  is 
the  assumption  that  certain  materials  are  capable  of  undergoing  small  plas- 
tic (permanent)  as  well  as  elastic  (recoverable)  strains  at  each  loading 
increment.  Mathematically,  the  total  strain  increment  is  assumed  to  be 
the  sum  of  the  elastic  and  plastic  strain  increments,  i.e.. 


de.  . = de.  , + de . , 
ij  ij  ij 


where 

de^j  = components  of  the  total  strain  increment  tensor 
E 

de^j  = components  of  the  elastic  strain  increment  tensor 
P 

de^j  = components  of  the  plastic  strain  increment  tensor 

6.  Within  the  elastic  range  the  behavior  of  the  material  can  be 
described  by  an  elastic  constitutive  relation  of  the  type 


de  = A (a  ) da,  . 
ij  ijkl  mn  kl 


where 


A. ,,  , = material  response  function 
ijkl 

dOj^^  = components  of  stress  increment  tensor 


The  behavior  of  the  material  in  the  plastic  range  can  be  described  within 

the  framework  of  the  generalized  incremental  theory  of  plasticity.  The 

1 P* 

mathematical  basis  of  the  theory  was  established  by  Drucker  by 


* Raised  numerals  refer  to  similarly  numbered  items  in  the  References 
at  the  end  of  the  main  text. 


A3 


introducing  the  concept  of  material  stability  with  the  following  impli- 
cations : 

(a)  yield  surface  (loeuding  function)  should  be  convex  in  stress 
space 

(b)  yield  surface  and  plastic  potential  should  coincide  (this  is 
referred  to  as  associated  flow  rule) 

(c)  work  "softening"  should  not  occur 

These  three  conditions  can  be  summarized  mathematically  by  the  following 
inequality 


do  de^  ^0  (a8) 

^0  1 J 

The  above  conditions  allow  considerable  flexibility  in  choosing  the  form 
of  the  loading  function  (j()  for  the  model,  which  serves  as  both  a yield 
surface  and  plastic  potential.  For  isotropic  materials  the  yield  surface 
may  be  expressed,  for  example,  as 


i (Ji 


(A9) 


where 


J,  = o = first  invariant  of  stress  tensor 
1 nn 


•^2  " 2 ^ij^ij 


second  invariant  of  stress  deviation  tensor 


S^j  = o^j  - (J^/3)  i5^j  = stress  deviation  tensor 


6.,  = 

ij 


= Kronecker  delta 


i = J 
(o  i j 
a hardening  parameter 


The  hardening  parameter  k , in  general,  can  be  taken  to  be  a function  of 


plastic  strain  tensor  e.. 

ij 


The  yield  surface  of  Equation  A9  may  expand 
or  contract  as  the  hardening  parameter  k increases  or  decreases,  respec- 
tively (Figure  Al). 

7.  Conditions  a,  b,  and  c above,  taken  in  conjunction  with 
Equation  A9j  result  in  the  following  plastic  flow  rule  for  isotropic 
materials 


Al+ 


if  ^ = 0 


de^  = 
ij 


if  1$  < 0 


where  dA  is  a positive  scalar  factor  of  proportionality,  which  is  non- 
zero only  when  plastic  deformations  occur,  and  is  dependent  on  the  partic- 
ular form  of  the  loading  function. 

Elastic  strain  increment  tensor 

8.  For  isotropic  elastic  materials  strain  increment  tensor 
(Equation  AT)  takes  the  following  form 


where 


elastic  bulk  modulus 


(All) 


G = elastic  shear  modulus 

The  bulk  and  shear  moduli  can  be  functions  of  the  invariants  of  stress 
tensor.  Accordingly,  we  will  assume  that  K = K (J^  ’ "^2  ’ 

G = G (J^  , Jp  , J^)  where  is  the  third  invariant  of  stress  deviation 

tensor.  Equation  All  can  be  written  in  terms  of  the  hydrostatic  and 
deviatoric  components  of  strain  and  stress  increment  tensors,  i.e.. 


^ 3K  (J^  , J^,  J3) 

de? . = dS . . 

2G  (J^  , , J3) 


(A12) 


(A13) 


de  = increment  of  elastic  volumetric  strain 
K-K. 


de.,  = elastic  strain  deviation  increment  tensor 


In  order  not  to  generate  energy  or  hysteresis  within  the  elastic  range, 
the  elastic  behavior  of  the  model  must  be  path  independent.  The  material 
should  then  possess  a positive  definite  elastic  internal  energy  function 
(W)  which  is  independent  of  stress  path.  The  strain  energy  function  can 
be  written  as 


W = / 


ij 


0 

a. 

f 

0 

J 


a.  . de. , 
ij  ij 


= (S  + i J 6 ) 

h ij  3 1 ij 


dJ, 


dS 


±- — 6 

9K  (J^  , Jg  , J3) 


ii. 


2G  (J^  , Jg 


S., 

= T IT dJ  + ^ ^ 

0 9K  \j  , J , JJ  ^ 0 2G  (J^  , , J3) 


= / 


d(J^)‘ 


dJ. 


^3^ 


S . . dS . . 

ij  iJ 


(AII4) 


0 18K  (J^  , , J3)  0 2G  (J^  , , J3) 

In  order  for  W to  be  independent  of  stress  path,  the  integrals  in 
Equation  AII4  have  to  depend  only  on  the  current  values  of  and 

Therefore,  the  bulk  and  shear  moduli  have  to  be  expressed  as 


K = K 
G = G iJ^) 


(AI5) 


Further,  K and  G must  always  be  positive.  Since  during  elastic  defor- 
mation the  hardening  parameter  (k)  is  constant,  the  bulk  and  shear  moduli 
can  also  be  expressed  as 


K = K (J^  , k) 
G = G (“2  , k) 


(AI6) 


A6 


Plastic  strain  increment  tensor 


9.  The  plastic  strain  increment  tensor  is  given  by  Equation  AlO 

where  the  loading  function  |$  is  given  by  Equation  A9»  The  hardening 

function  in  Equation  A9  could  be  taken  as  being  equal  to  plastic  volu- 
P 

metric  strain  e,  , , thus 

kk  ’ 


K = 


(A17) 


The  use  of  Equation  AIT  will  allow  the  yield  surface  to  expand  as  well  as 
to  contract.  Figure  Al.  The  plastic  loading  criteria  for  the  function 
1$  are  given  as 


^ > 0 loading 
9 

dX  = . ( =0  neutral  loading  (Al8) 

30ij 

^ < 0 unloading 


Because  de . . =0  during  unloading  or  neutral  loading,  as  well  as  for 
< 0 , Equations  All  through  A13  are  used  to  determine  the  purely  elas- 
tic strain  changes.  The  prescription  that  neutral  loading  produces  no 
plastic  strain  is  called  the  continuity  condition.  Its  satisfaction  leads 
to  concidence  of  elastic  and  plastic  constitutive  laws  during  neutral 

n 12,13 

loading.  ’ 


10.  Like  the  elastic  behavior,  the  plastic  stress-strain  relation 
can  be  expressed  in  terms  of  the  hydrostatic  and  deviatoric  components 
of  strain.  Applying  the  chain  rule  of  differentiation  to  the  right-hand 
side  of  Equation  AlO  yields 


AT 


or 


/ 


de 


dc". . d),r2jL 


<5.  , + 


1 3 


ij  3J,  ij  f=-  [=■  ij 


(A19) 


Multiplying  both  sides  of  Equation  A19  by  6..  gives 


(A20) 


The  deviatoric  component  of  the  plastic  strain  increment  tensor  (de  ) 

^ tj 


can  be  wri hten  as 


P PIP 

de . . = de  - — de  6 . . 
ij  ij  3 kk  ij 


Substitution  of  Equations  A19  and  A20  into  Equation  A21  yields 

J’  = dX  3j( 


(A21) 


de. 


ij 


(A22) 


11.  In  order  to  use  Equations  A20  and  A22,  or  Equation  A19,  the 
proportionality  factor  dX  must  be  determined.  This  can  be  accomplished 
in  the  lollowing  manner.  From  Equations  A9  and  A17  the  total  derivative 
of  becomes 


l£ 


1 " 2-|/^3-#2 


S . , dS . . + de^.  = 0 


ij  IJ 


3e 


kk 


(A23) 


kk 


In  view  of  Equations  A12,  A13,  and  A20,  Equation  A23  becomes 


(A2it) 


Substituting  Equation  a6  into  Equation  A2'k  results  in 


® '""kk  ' IjT  ^ 7§-  <^»ij 

V'>2 


Ij  ^rzr  Ij 


-3  dA  (A25) 

kk 


3"  <»Sk  If  * — ^ ' 3K  dJ,  If  * ^ ii-  S.  dk^ , 

iJ,  if,  f 

^ ^^kk 


(A26) 


P P 

Substituting  the  values  of  and  de^^  from  Equations  A20  and  A22, 

respectively,  into  Equation  A26,  we  obtain 


3K  de,.  If  + ^ S,.  de,  = 9K  dA  (|f )"  + G dA 


'kk  dJ  ^ S.jde.j  = 9KdA  ( : 

" K ^¥2  " 


^ ^"kk 


Solving  for  dA  , we  obtain 


(A27) 


3K  If  de  + — — S.  de 

kk  r=^  r=:r  ij  ij 


lfj2  sYf 


« <lf  ’ If  !V 


(A28) 


Total  strain  increment  tensor 

12.  The  total  strain  increment  tensor  can  be  obtained  by  combining 
Equations  A6,  All,  AI9,  and  A28;  thus, 


dJ  dS . 

^ 9k"  '^ij  2G 


3K 


de  + S de 


9K  (|^)  + G 


-li  T _ -3  3|^  ii 


'■^1  a/  ^ 

kk 


fliL  5 + _1 3jLs 

3Jn  ij  !=-  r=-  i 


(A29) 

J 


Similarly,  the  stress  increment  tensor  can  be  written  as 


da  = K de,  , 6.  ^ + 2G  de. , 
IJ  kk  ij  ij 


3K  de  + — S de 


9K  (|f-)  + G 

3Jl 


,3j_  li_ 


ay  j 


a p 

1 ^"kk 


(A30) 


iLfiLl 
)iT  2i[j~-> 


Equation  A29,  or  Equation  A30,  is  the  general  constitutive  equation  for 
an  elastic-plastic  isotropic  material.  To  use  these  equations  it  is  only 
necessary  to  specify  the  functional  forms  of  K , G , and  ^ . 
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